Application 



For 



United States Letters Patent 



SPECIFICATION 



TO WHOM IT MAY CONCERN:- 



BE IT KNOWN, THAT I, Robert G. Komarechka, residing at 545 Granite 
Street, Sudbury, Ontario, P3C 2P4, Canada, a citizen of Canada respectively, have invented 
or discovered certain new and useful improvements in:- 



METHOD OF DISPLAYING 
THREE-DIMENSIONAL VECTOR ORIENTATIONS 
ON A TWO-DIMENSIONAL SURFACE 



of which the following is a specification. 



TITLE OF THE INVENTION 



METHOD OF DISPLAYING 
THREE-DIMENSIONAL VECTOR ORHSNTATIONS 
ON A TWO-DIMENSIONAL SURFACE 

FIELD OF THE INVENTION 

The present invention relates, in general, to data processing and color-coded 
raster rendering of information, and more specifically, to a method for displaying the 
orientation of three-dimensional vector data on a two dimensional display or surface. This 
methodology has specific applications in the field of geophysical surveys. 

BACKGROUND OF THE INVENTION 

The term "ferromagnetic" used hereinafter shall constitute the inclusion of 
antiferromagnetic and ferromagnetic materials. 

The term "Currie temperature" or "Currie Point" used hereinafter shall also 
encompass the Ne'el temperature as would be used in the context of antiferromagnetic and 
ferromagnetic materials. 

In the study of vectors within vector fields it is often useful to correlate 
information that relates to the magnitude, orientation and spatial relationship of the 
individual vectors within said field. 
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While numerous effective methodologies exist for displaying the magnitude 
of individual vectors, the ability to observe commonality of individual vector orientation 
within vector fields containing thousands of individual vectors on a surface or within a 
volume can be difficult to discern. 
5 The applicant is aware of several attempts in prior art, which relate 

particularly to the visualisation of three-dimensional vector fields. Reference may be had by 
turning to the following publications hereinafter further referred to in superscnpt ; 

Reference 1 : "Visualizing the Magnitude and Direction of Three-Dimensional 
Vector Fields" by Timothy J. Peter in AP-S International Symposium Digest, IEEE 
1 0 Antennas and Propagation Society, Volume 3, 1 993, pages 1 1 1 8 to 1 1 2 1 . 

Reference 2 : "Visualization of Three-Dimensional Vector Fields and 
Functions" by J. Simkin, IEEE on Visual 3D Fields Digest, 1995/205, pages 1/1-2; and 

Reference 3 : "Vector Field Animation with Maps" by Yamrom, IEEE 
Computer Graphics and Applications, Volume 15, Part 2, pages 22 to 24, March, 1995; 
15 Reference 4 : "Combined 3D Interpretation of Airborne, Surface and Borehole 

Vector Magnetics" of Mueller, E. L., Morris, W. A., Killeen, P.G., and Balch, S. at the 
McConnell Nickel Deposit, Proceedings of Exploration 97: Fourth Decennial International 
Conference on Mineral Exploration, edited by A. G. Gubins, 1997, p. 657-666. 

Reference 5 : "Mapping Remnant Magnetism Using the Local Phase" of 
20 Thruston, J. B., in Geophysics, Volume 66, No. 4, July-August, 2001, p. 1082-1089. 

Current methodologies for indicating vector orientation consist of 
streamlines, hedgehogs, stream-ribbons, stream-polygons and textures maps 3 . In all these 
methods there are difficulties in observing the vector orientation orthogonal to the 
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observation surface or in selection of unique desired orientations. Figure 5 of the drawing 
sheets illustrates a vector grid that shows this effect A novel texture mapping procedure is 
hereby disclosed that overcomes many of the shortfalls of the earlier methods. 

Many applications exist for the use of an effective display of vector 
5 orientation including studies in fluidics, electromagnetic fields, finite element analysis, 
magneto hydrodynamics, plasma physics, astrophysics magnetic and digital elevation 
models (DEMs). We will present an embodiment in which the application of the present 
method is illustrated in the field of geophysical geomagnetic mineral exploration. It is to be 
recognized that this embodiment is but only one example of the usefulness of this Vector 

1 0 Orientation Visualization Method. 

Geomagnetic surveys have many uses and are of particular interest in mineral 
exploration, as they may identify areas likely to contain some types of valuable deposits. 
They can also identify various lithologies, and provide some types of structural information. 

A common method of performing such a survey is to measure the magnetic 

15 field at various points along the Earth's surface while recording the magnitude variations in 
the magnetic field strength. Although dip needle surveys have been used in the past, until 
recently, the orientation of the magnetic field was not generally measured. The magnitude 
of variations in the readings was primarily due to variations in the magnetic susceptibility of 
the constituent minerals within the rock, and to a much lesser degree, the variation of the 

20 orientation of the remnant field. Such surveys, while giving useful data, which related 
primarily to lithological effects, disregarded the information of the orientation of the 
magnetic field. 

Recently, due to the introduction of tri-axial gradient magnetometer surveys, 
it has become possible to retrieve vector oriented gradient data. In addition, tri-axial 



fluxgate total magnetometer surveys can also be used to measure the orientation of the total 
magnetic field vector. Tri-axial fluxgate magnetometers have also been used in along a 
borehole to effectively delineate the partial outline of a dipole about a ferromagnetic ore 
body 4 . 

The ability to effectively display the orientation of the magnetic vectors over 
a surface or within a block of rock as by borehole surveys can greatly add to the geologic 
understanding of the survey area and can lead to a better understanding of the location and 
extent of several varieties of ore deposits. 

To understand the benefit of knowing the orientation of the vectors within the 
survey area it is important to understand the means by which variations in the orientation of 
the vectors can occur. After a high temperature thermal event, certain minerals in the earth 
will heat enough beyond their Currie Point Ne'el temperature to lose their magnetism. 
When the temperature lowers, these minerals regain their magnetism and will now have a 
thermal remnant magnetic field aligned with that of the earth at the time of going through 
their Currie temperature. Since the earth's magnetic field is known to change directions and 
flip at various times, the differently oriented remnant magnetic vectors of these 
remagnetized minerals can indicate spatially segregated intrusives of contemporaneous age, 
multiple pulsed intrusives within a similar rock type, the extent of their thermal aureoles and 
the presence of folds and faults. Similarly alteration due to chemical -known as chemical 
remnant magnetism, CRM- and pressure changes -known as piezoremnant magnetism- can 
result in the formation of new minerals, which will retain a remnant orientation of the 
earth's field orientation at that time -known as detrital remnant magnetism (DRM) caused 
by the quiescent sedimentation of fine magnetic material aligning with the earth's magnetic 
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field at the time of sedimentation.- DRM is another remnance whose variation can be used 
to recognize later faults and folds. 

For the purpose of this embodiment a magnetic survey can be performed 
using either a tri-axial fluxgate magnetometer or tri-axial gradient magnetometer. The 
5 vertical gradient can also be calculated if both horizontal gradients and the total magnetic 
field are known. These types of surveys produce the magnitude of the magnetic field for 
each of the three-orthogonal directions. The data is collected for process and analysis over 
the survey area by establishing a two-dimensional grid over the area of interest. 
Measurements of the magnetic field are taken at specific points on the grid either by aircraft 

10 or by an operator on the ground working along a cut or demarked grid. . Each collection 
point is identified as a point on the two-dimensional grid, and will have an associated three- 
dimensional vector representing the magnetic field at that point. Alternatively, if the data 
should be collected along a borehole as a function of its distance from the collar, the location 
of the three-dimensional vector can be represented from the surveyed borehole information 

15 as a point in 3D space with x, y, z co-ordinates. The x,y co-ordinates can match the two 
dimensional grid of the surveyed area of interest mentioned earlier. 

Peters Reference 1 : This reference describes a method of representing an 
electromagnetic vector field defined over a three-dimensional surface. A triangular grid 
represents the surface, and the method uses color shaded magnitude contours and directed 

20 vector lines to characterize the vector field. Referring to the Figures, it may be seen that in 
this reference, the direction of a vector field is represented using arrows, and the magnitude 
of the vector field is displayed using colors. Such visualization is suited to vector fields for 
which the direction changes gradually and continuously. The magnitude of the vector field 
can then be visually identified by the intensity of the colours used. 



Simkin Reference : This reference provides an overview of three- 
dimensional visualization, and only contemplates displaying vectors using arrows. The 
direction of the vector is shown using arrows while the length of the arrow visualizes the 
magnitude, In other words, the longer the arrow, and the higher the magnitude of the vector. 

Yamrom Reference 3 : This reference describes a method of animating vector 
fields and discusses the subject matter of displaying a 3D-vector field by changing the 
texture of the viewed surface with respect to time. 

SUMMARY OF THE INVENTION 

According to an aspect of the present invention, there is provided a method of 
displaying three-dimensional vector field information on a two-dimensional display* After 
collecting the information from an area of interest, the Cartesian vectors are transformed into 
spherical co-ordinates. A colour model such as Munsell, HSV, HLS, or YIQ, is modified 
into a spherical colour model. In the embodiment of the modified HSV model, the spherical 
model is developed such that Hue is distributed evenly along the equator while the Value 
varies from white at the North Pole to black at the South Pole. Each of the vectors are 
plotted emanating from the origin of the coloured sphere using their spherical co-ordinates. 
An HSV colour is assigned to the area of the vector on the map based on the colour 
occurring at the corresponding piercing point of the vector with the designated HSV 
spherical colour model surface. This HSV colour is then used to fill the area or volume of 
influence of the vector. In the case of a monitor display the HSV values are converted to 
their corresponding RGB values. Ability exists using this model to define a particular 
orientation range of azimuth and dip such that the whole colour gamut may be displayed 
within a narrow range and hence further delineate variations in orientation. 



The utility of this methodology includes but is not limited to geophysical 

surveys. 

For purposes of clarification and to provide a full disclosure, it was deemed 
absolute that color drawing be included in the present application in order to fully 
understand the scope of inventiveness. 



BRIEF DESCRIPTION OF THE DRAWINGS 



These and other advantages of the invention will become apparent upon 
reading the following detailed description and upon referring to the drawings in which:- 

Figure 1 is a flowchart outlining an embodiment of a method of data 
collection and display of three-dimensional information using a two-dimensional display. 



Figure la is a flowchart outlining an embodiment of a method of data 
collection and display of three-dimensional information using a two-dimensional display 
and indicating the variables used to describe the data and its manipulation. 

Figure 2 is a schematic diagram of a grid of irregular outline established to 
define the collection sites of the three-dimensional information. 

Figure 3 is an example of a standard base station calibration. 

Figure 4 is an illustration of the definition of a spherical co-ordinate system 
based on Mathematical Conventions. 



Figure 4a is an illustration of the definition of a spherical co-ordinate system 
based on Geological Conventions. 



Figure 5 is a schematic diagram of a two-dimensional grid surface defined for 
locating the collection sites of the three-dimensional information (3D vectors). 
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Figure 6 illustrates the notation used in the orientation of the axis and co- 
ordinate systems. 

Figure 7 shows an HSV colour coded sphere sectored in planes of equal area, 
each denoting a specific HSV colour with the 3D vectors being oriented so as to emanate 
from the origin to intersect a uniquely coloured surface of the sphere. It should be noted that 
this sectored sphere is presented in this illustration for the basis of clarity of methodology. 
In reality a continuum of millions of colours will be more frequently used. 

Figure 8 is a table illustrates azimuth and dip angles of the 3D vectors as per 
the grid in Figure 2. 

Figure 9 is a table illustrating assigned HSV colour values to the 3D vectors 
as per the grid in Figure 2. 

Figure 10 illustrates a further endowment of this methodology that allows for 
the ability to illustrate a scattergram distribution of vector orientation on a projection. 

Figure 10a illustrates a section taken from Figure 10 further expanding the 
details of the orientations of the vectors. 

Figures 11, 11a and lib illustrate an endowment of this methodology that 
shows how a range of vector orientations may be sectored out and given higher visual 
discrimination to further emphasis subtle changes in vector orientation. 
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Figure 12 illustrates the described methodology of colour coded vector 
orientation as being applied to voxels located at sampled points along the 3D survey 
trajectory of a curvilinear borehole. 

Figure 13 HSV colour assignment illustrates the relationship of the observed 
5 geological spherical co-ordinates and the development of the modified HSV spherical colour 
model. 

While the invention will be described in conjunction with illustrated 
embodiments, it will be understood that it is not intended to limit the invention to such 
embodiments. On the contrary, it is intended to cover all alternatives, modifications and 
10 equivalents as may be included within the spirit and scope of the invention as defined by the 
appended claims. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

In the following description, similar features in the drawings have been given 

similar reference numerals. 

Turning to Figure 1, a flowchart outlining an embodiment of a method of 

5 displaying three-dimensional information, such as three-dimensional vectors, using a two- 
dimensional display, is shown. More specifically, the present embodiment is directed at the 
collection of three-dimensional magnetic field data 14 in an area of interest 10 on the surface 
of the Earth from a plurality of mobile stations to extract the residual magnetic field 
orientation of the underlying rocks. Said data is then calibrated 16 at a base station 12, 

10 converted to spherical co-ordinates 18, assigned colour values 20 then plotted 22 to a two- 
dimensional surface. A further detail of this flowchart and the variable names used to 
describe this data is shown in Figure la. Although the present embodiment is directed at a 
method of displaying three-dimensional magnetic field vectors for use in the study of 
geophysics, thermodynamics, fluidics, plasma physics, magneto hydrodynamics and the 

15 like, the method may also be used to display vectors representing any three-dimensional 
information in a two-dimensional display. 

Initially, in the case of a surface ground based geomagnetic survey the area of 
interest is selected and defined by setting up a grid over the area form which the magnetic 
field data is to be collected 10, in the case of an aircraft survey, location of the flight lines 

20 would be obtained by utilizing rapid precision Global Positioning System (GPS) equipment. 
Figure 2 illustrates an example of a grid 100 including a carefully laid out baseline 102 
which is located centrally within the grid and more or less parallel to the feature desiring 
investigation and dividing the grid area into two parts 104 and 106. This baseline 102 is cut 
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free of brush and trees so as to allow for a line of sight to measure out the locations of the 
lines 112. Due to the trend of the baseline in Figure 2 being in a more or less South-North 
direction, each side of the baseline is segregated in an Eastern and Western section 
respectively. Lines 112 perpendicular to the baseline 102 further define the grid 100 and are 
demarked from the baseline via surveying or by the use of a turning board. Said lines 112 
are then blazed or line cut and picketed at regular station intervals for later measurement 
with the instrument. In the exemplary illustration presented, the lines 112 are shown 100- 
meters apart with measuring stations 116 at 25-meter intervals along the line emanating 
from the baseline 102. For computational purposes, an outline of rectangular area 114 is 
defined around each of the measuring stations 116 as follows: In the direction of the lines 
112, the opposite sides of the rectangle are located at the midpoint of each of the consecutive 
measuring stations 116, which in the example case would be 12.5 meters. In our example, 
12.5 meters would be used to define the on-line rectangle dimension of the measuring 
stations at the end of each of the lines 112. Orthogonal to this, the distance between sides of 
the rectangle can be arbitrarily set between the previously determined online side-dimension 
as a minimum and the distance between lines 112 as a maximum Various algorithms are 
known in the art for calculation of values between lines 112. 

At each of the measuring stations 116, the magnetic field data, the time of the 
measurement and the station location is recorded within the grid 100. 

In order to distinguish measurements from the various measuring stations 
116, the measuring stations are denoted as LN(orientation)+(distance)(section). This is a 
standard notation in the exploration industry in North America and is further explained 
below. The 'LN' values represents the measuring station's line location with respect to a 
selected origin of the baseline 102 with L0+00 representing the measuring station at the 
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origin on L00+00. In our example, since the lines 112 are orthogonal to the baseline and 
100-meters apart, n-BLdistance/100, 'distance' being the distance along the baseline 102 
from the selected origin. As one progresses in a northward 'orientation', the lines 112 are 
sequentially labelled LOO, LOIN, L02N .... LN with (N+l) being the number of lines 112. 
The 'distance' value represents the distance the measuring station 116 is located away from 
the baseline 102 and the 'section' value represents the section 104 or 106 in which the 
measuring station is located. Therefore, for a measuring station in the Eastern section 104 
located 50-meters northward away from the baseline 102 on the perpendicular line 
corresponding to L04N, this measuring station may be denoted as L04N+50E. In the grid of 
Figure 2, the number of baseline measuring stations 116 is indicated by the tick marks along 
the lines 112. Further reference to Figure 2 illustrates the designation of the measuring 
stations 116 along the baseline. In respect to those measuring stations 116 occurring along 
the baseline the notation BL(orientation)+(BLdistance)(direction) is used. 

After the grid 100 has been defined, a magnetically quiet area, close to the 
grid, is selected so that a tri-axial magnetometer base station 118 may be located (step 12) 
and appropriately positioned to collect signals which are used to assist in the calibration of 
the signals collected at the measuring stations 116. The base station 118 instrument's clock 
is set to match the clock of the portable mobile tri-axial instrument. This will be described 
in more detail below. 

Prior to collecting the data (step 14) the operator initially records, from a 
distance, x,y,z readings in a magnetically quiet area with the tri-axial magnetometer, this 
could be a tri-axial fluxgate instrument suspended at the proper height and orientation by a 
remote nonmagnetic device. A second set of readings is then taken in the same location 
with the operator holding the instrument in the normal manner to take a reading. These two 
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readings are recorded so that any magnetic effect from the presence of the operator can be 
removed. This procedure may be undertaken to confirm a consistent operator null 
correction 

Having undertaken the above operator nulling procedure, the operator walks 
with the instrument along the baseline and lines such that the individual bypasses each 
measuring station at least once. At each measuring station, the individual records and stores 
magnetic field data in the form of magnetic field vector information along with the station 
location and the time the measurement is recorded. Care is taken in the positioning and 
orientation of the instrument, aided by a sighting and levelling device to maintain a 
consistent x,y,z orientation of the probes. The magnetic field vector information collected 
includes information corresponding to the magnitude and direction of the magnetic field at 
the measuring station location. This may be achieved by simply writing down the magnetic 
field data on a notepad or by entering such information into a data collector for later 
dumping on a computer. The recorded signals Xj,Yi and Zi may be represented as X m , Y m 
and Zm with " m " representing the measuring station designation. 

At the end of each days work after the individual has collected the magnetic 
field data from each of the measuring stations visited, the individual retrieves calibration 
information from the base station corresponding to the time period over which the magnetic 
field data was being recorded. After retrieving the calibration information from the base 
station, the calibration information and the recorded signals are stored in a computer and the 
recorded signals are calibrated. This is a standard procedure known to those skilled in the 
art. 

The calibration process allows for the magnetic field data retrieved from the 
measuring stations to be a truer indication of their actual magnitude and direction. It will be 
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understood that the magnetic field data may be affected by external conditions such as solar 
storms or diurnal variations. 

Turning to Figure 3, a sample display of calibration information is shown. 
Base station averages, Xbavg, Ybavg, and Zba V g, are calculated by averaging the magnitude of 
the stable magnetic field vector, in the X, Y and Z directions at the base station over the time 
period of the survey. Since the base station is assumed to be located in a magnetically quiet 
location away from vehicular traffic, power lines, transmitter sites, etc. it may be assumed 
that any recorded deviations from the base station stable average indicates that any 
measurements taken at that instance in time have been affected by the same external 
conditioa For instance, at time = 6:15, as indicated on Figure 3, the X value measured at 
the base station is seen as a spurious X b positive peak above the pre-calculated average, X b 
and therefore a base station calibration correction value X teses tation equals the absolute value 
between X b and the pre-calculated average, Xba Vg . 

After selecting a first-recorded signal X m , Y m , 2^ from the mobile station, the 
time at which the signal was recorded is retrieved. The calibration values for this time is 
then calculated as explained above. If the reading was taken at 6:15, the base-station 
calibration values Xhasestation6:i5, Ybasestation6:i5, and Zbasestation6:i5 are retrieved and subtracted 
from the measuring station values X m , Y m and Zm corresponding to the first recorded signal 
in order to perform a first step of calibrating the measuring station value. After the first 
reading is calibrated, the magnetic field data from a second reading is selected and calibrated 
by retrieving the time measurement from the magnetic field data and calculating and 
retrieving the calibration values corresponding to that time measurement. This step is 
repeated until all of the magnetic field data from each of the measuring stations is calibrated. 
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Therefore Xcaiibrated equals X m - XbasestationCtime) which is equal to X m - |Xbav g - Xb(time)l- Values 
for Yetted and Zested are calculated in a similar manner. 

After the values Xcaiibrated, Ycaiibrated and ZcaHbrated for each measuring station 
have been calculated. A second step may be performed to remove the geomagnetic field 
vector from these values in order to obtain a residual magnetic value from the magnetic field 
data. 

This is achieved by calculating the Earth's magnetic field using the standard 
International Geomagnetic Reference Field (IGRF) procedures to determine the theoretical 
values with respect to the area of interest at the time of the survey These theoretical values 
obtained may be represented as Yigrf and Zjgrf. 

The residual magnetic data values Xresiduai, Yresiduai and Zresiduai may then be 

Calculated from the equations Xcaiibrated" Xigrf, YcaHbrated - Yigrf, and Zcalibrated - Zigrf, 

respectively. 

After calculating the residual magnetic values, further methods known to the 
trade 5 may be performed which causes the induced magnetic field value to be removed from 
the residual magnetic data values to obtain remnant field values at the measuring stations. 

In order to calculate the remnant values Xremnant, Yremnam and Zremnam, the 
equations Xresiduai- Xinduced, Yresiduai - Yinduced, and Zresiduai - Z^ced, are respectively 
determined. The remnant value represent the magnetic field data values were frequently 
acquired at the time the ferromagnetic rock forming minerals were formed or went through a 
Currie temperature episode of its ferromagnetic minerals. This remnant vector data can 
allow for the discrimination as outlined earlier in the claims 79(c). After calculating either 
the residual or remnant data these x, y, z Cartesian values are transformed to mathematical 
spherical co-ordinates using the following equations: 
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Tmath = squareroot(X 2 + Y 2 + Z 2 ) where r is the radius of the unit sphere used 
to determine spherical co-ordinates; 

©(theta^th = arctan (Y/ X) where theta represents the azimuth angle, which 
is measured from the positive x-axis toward the positive y-axis theta and is a value between 
0° and 360°; 

0(phi)math = arctan ((squareroot(X 2 + Y 2 )/ Z) where phi is the angle measured 
from the vertical positive z axis counterclockwise toward the x,y plane and is a value 
between 0°to 180°. 

These mathematical spherical co-ordinates are illustrated clearly in Figure 4. 

The mathematical spherical co-ordinates of the magnetic vector at each 
measuring station may then be represented by (rmath, ©math, ^math). 

A conversion of the mathematical spherical co-ordinates (r^th, ©math, <J>math) is 
then translated into the geological co-ordinates (r g , ©g, Og) by the following equation: 
(r g , © g , 0>g) = [(rmath), (©math - 270), (90- Omam)] as clearly illustrated in Figure 4a. These 
values are tabled in Figure 8. 

After calculating the geological spherical co-ordinates for each measuring 
station, a unique colour based legend may be determined by assigning colour or patterns to 
the spherical co-ordinates. Several colour models or notations exist which have direct 
symmetry with respect to spherical co-ordinates such as red, green, blue (RGB); 
Commission Internationale de TEclairage (CIE), hue, saturation, value (HSV); hue, 
saturation, luminance (HSL), CIE XYZ, YIQ, Munsell, TekHVC, CIE LUV etc... The last 
three of these models allow for a more continuous perception of colour to the human eye. 
Although any of these colour models may be used in the preparation of a suitable colour 
model via the proper mathematical translation. In this embodiment a modified spherical 
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HSV colour model is used, where hue variation is assigned the azimuth angle. Axial 
symmetry about dipole axis may also be illustrated by having hue vary as dip. 

One example of assigning spherical co-ordinates to an HSV model may be 
realized by the following two equations, the use of each being dependent on the value of phi: 
When O>0 [H, S, V] « [(0g), (2 O g /7i), (1)], when <D<0 [H, S, V] = [(0 g) , (1), 
(2-2 Og/rc)], Figure 6 illustrates this colour assignment. 

The above assignment relates a unique HSV colour value to the direction of a 
three dimensional vector corresponding to measured and calculated values from the 
measuring stations which can be applied to the area of interest. 

Figure 9 graphically illustrates this principle. The particular range of values 
for the HSV colour values defines the limits of the above mathematical expressions. 
However, other parametric vector assignments can be used for different colour models or for 
the creation of a new unique colour model or even textural patterns. 

After each of the spherical co-ordinate vectors have been associated with an 
HSV value, see the table of Figure 10, the colour representing the HSV value is then applied 
as a fill colour to a rectangular area centred at each of the corresponding measuring station 
locations along the lines with respect to the overall grid. Figure 5 shows the vectors as 
arrows occurring at the centre of a rectangular cell around the location of their measurement 
point illustrating the difficulty in detennining change in magnitude, verses a change in 
foreshortened length due to a change in dip. 

Figures 8 and 9 provide two-dimensional displays of the magnetic field data 
in accordance with an embodiment of the present method. For enhanced clarity, in Figures 5, 
8 and 9, the lines of the grid are not shown but lie along the centre of each grid cell. 



20 



These coloured areas located about each of the measuring stations may also 
be input into any of the many standard contour algorithms that exist for plotting of these 
coloured pixelated areas to produce a colour map with gradational interstitial colours applied 
between the lines and covering the whole area. These interstitial colours could represent, 
along observed trends, the shortest distance between the two colours as illustrated on the 
spherical colour model. 

With the advent of computer processing capabilities, the unique display of 
this vector orientation data lends itself to a number of useful statistical presentations. In 
Figure 10 we see that an area of the colour coded data grid area can be selected (Figure 10a) 
and the orientation of all vectors within this area can be plotted on an orthogonal grid as a 
function of dip angle versus Azimuth (Figure 10b). This can allow for the fine cognitive 
discrimination of very subtle orientation regimes. 

Figure 11 shows another attribute of this vector orientation presentation, 
which allows for the enhanced discretion of specific selected orientation. In this scenario a 
specific range of dip and azimuth angles are selected. This selection is illustrated in the 
highlighted quadrant of the sphere (Figure 11a) or the rectangular area expressed on the 
orthogonal dip verses azimuth plot (Figure lib). A new colour chart is then created that 
contains the full complement of colours as the original model (Figure 1 1c). The appropriate 
cells that contain the dip and azimuth values within this range are displayed on the plan view 
of the grid area utilizing the fiill gamut of colours available. 

Additionally the magnitude of the vectors may also be shown on the coloured 
illustration either as dashed black and white line contours or as a Digital Elevation Model. 

Alternatively to the ground acquisition method described earlier, an aircraft 
flying over the grid may retrieve the above information. After the grid has been determined, 
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the aircraft flies over a prescribed calculated grid lines to record the magnetic vector 
information at regular time intervals corresponding to a specified distance given the airspeed 
flown. In order to ensure that the aircraft follows the calculated grid lines, a GPS is be used 
to guide the aircraft. In this aerial data collection no baseline or ground demarcation of lines 
5 is required. 

The magnetic field vector information is collected via instruments on the 
aircraft as they fly over the measuring station and stored on record keeping equipment, such 
as a data logger, located on the aircraft. Calibration and corrections to the data are carried 
out later after each day's flying. 
10 An example of the 3D versatility of this colour coded vector orientation 

methodology is also given in Figure 12 where the method illustrates vector orientation 
sample points collected along a borehole within a given volume by the use of coloured 
voxels. 
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